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C^ I numerically analyze Joule heating in large Bi2212 mesas while taking into account typical 
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thermal conductivities and their temperature dependencies of all the materials involved in the heat 
dissipation and its removal. Such mesas are used in experiments on THz-range radiation. The 
analysis shows that the temperature increases with bias current and is distributed unevenly along 
the mesas. The temperature of the mesa's middle part can even exceed Tc at sufficiently high 
bias. The non-uniform temperature distribution can possibly be important for synchronization of 



c/3 emission from different junctions in the mesas. The overall cur rent- volt age characteristics are also 

Cd calculated self-consistently showing a negative differential conductance in a wide range of currents. 
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T^ PACS numbers: 74.25.Sv, 74.50.+r, 74.72.Hs 

^ Keywords: Intrinsic Josephson junctions, THz radiation 
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I. INTRODUCTION 

The ac Josephson effect provides the basis for superconducting high-frequency detectors 
and tunable radiation sources with a voltage-to-frequency conversion parameter involving 
only fundamental constants. Josephson junctions made of high-temperature superconduc- 
tors (HTS) allow for such devices operating in a wide temperature range ~ 4 — 90 K and high 
frequencies including the THz band. There are not many continuous-wave tunable sources 
available at the low end of this band: 0.3 — 1 THz (see a review article, Ref. pj). The HTS 
Josephson junctions would be almost ideal radiation sources at these frequencies if they 
only had sufficiently high output power. Thin films of YBa2Cu307+5 (YBCO) have been 
the material of choice for making HTS Josephson junctions. The Fiske- and external-cavity 
resonances were clearly seen in such junctions that witnessed for the presence of Josephson 
radiation [21 13]. However, the radiation from the individual junctions is weak, usually a few 
tens of pW, which limits their utility for practical applications. 

Synchronization (or phase-locking) of many Josephson junctions is needed to boost the 
radiation power. A few YBCO Josephson junctions made on the same chip were demon- 
strated to phase-lock and thereby make the emission stronger [U |5]. The synchronization 
looks easier to achieve for the so called intrinsic Josephson junctions (IJ J's) which are densely 
packed inside single crystals of another popular HTS, Bi2Sr2CaCu208+5 (Bi2212) [6|. The 
high-frequency radiation from these junctions was seen in several works [7, 8J. However, the 
radiation was not strong either which indicated an absence of junctions phase-locking. 

A relatively intense THz-range emission has only recently been detected from large-area 
and high mesas etched into Bi2212 single crystals [9]. A number of models have been pro- 
posed to explain the phase locking of I J J's in these mesas [T0Hl2]. Yet stronger radiation has 
been measured pi3j at a current much higher than that in the original work of Ref. [9]. This 
current corresponds to a back-bending region of the overall current- voltage (I-V) character- 
istics, i.e. the region with a negative differential conductance. Joule heating is apparently 
high at this bias and therefore requires analysis of the heat balance and cooling efficiency. 

In this work, I calculate the temperature distribution in large-area mesas. The cal- 
culations include realistic thermal conductivities of all the materials involved in the heat 
generation and removal, as well as their temperature dependencies and anisotropy. In other 



words, I numerically solve the non-linear diffusion equation 

V(-A;(T)VT) = p(T)j2 (1) 

for mesas of typical geometry. Here, k{T), j, and p are the anisotropic thermal conductiv- 
ity,the current density, and the anisotropic resistivity, respectively. 

The calculations show that the temperature distribution in the large-area mesas is non- 
uniform for a wide range of currents. The non-uniformity increases with current. At suffi- 
ciently high bias, the temperature in the center even exceeds the superconducting critical 
one, Tc, while the rest of the mesa is still superconducting. The appearance of the normal- 
state regions can be important for synchronization of radiation from many IJJ's at high 
bias [13j. 

II. GEOMETRY OF MODEL MESA 

Fig. [T] schematically shows a mesa on top of a Bi2212 single crystal. The crystal is taken 
to be 1 X 1 mm^ large and either 20 or 40 /im thick. The single crystal is attached to a 
sapphire substrate with a 20-/im-thick polymethylmethacrylate (PMMA) glue layer. Due 
to the very high thermal conductivity of sapphire, the substrate surface is assumed to be at 
the bath temperature Tq. There is no cooling gas around the sample meaning that the heat 
dissipated in the mesa can only escape through the underlying single crystal and the PMMA 
layer. The calculations are made for one and the same length of the mesa / = 300 ;um. The 
mesa width w (the height h) is either 50 or 100 /im (1 or 2 yum). A 100-nm gold thin film 
is deposited on top of the mesa and the front part of the single crystal making up the bias 
electrodes (see Fig. [I]). There is no interim resistive layer between the electrodes and Bi2212. 
All the simulations have been made using the finite-element software package COMSOL [13]. 

III. MATERIALS PARAMETERS 

For calculations, I use common materials parameters found elsewhere. The anisotropic 
thermal conductivity k{T) of Bi2212 single crystals was measured in a number of works, 
predominantly along the a6-plane [T3HIH]- To the best of my knowledge, there exist only two 
measurements of the thermal conductivity in the c-axis direction [T^l ED] • It was found that 




FIG. 1: (Color online) Schematic view of a mesa used in the calculations. The height /i = 1 or 2 
/im; the width it) = 50 or 100 /xm. The length I = 300 //m, the single-crystal thickness d and the 
glue thickness t are both 20 /im. Not in scale. 



the ratio kab/kc ~ 6 in a wide temperature range and increases to 8 at low temperatures [I9] . 
In the calculations, I use kab for pure Bi2212 adopted from Ref. [18] and simply assume that 
kc = kab/lO. 

The a6-plane resistivity of Bi2212 is taken to be linear in temperature: pab = 6{T — 
Tc,6T)T/75 [fiQm], where 6{x,6x) is the Heaviside step function smeared over interval 6x 
and Tc = 86 K is the superconducting critical temperature of the optimally doped Bi2212. 
The Heaviside function is used to mimic the transition to the normal-metal state when 
the local temperature exceeds Tc. The current-induced transitions within individual CuO 
planes [21] are ignored in these calculations. 

The c-axis resistivity pc below Tc can be obtained from the sub-gap resistance Rsg of 
IJJ's at different temperatures [231 [21] • Rsg is estimated in the limit /, V — )■ and is there- 
fore not affected by Joule heating. Rgg shows a tendency to saturate at low temperatures 
possibly reflecting the presence of weakly-conducting channels short-circuiting the interlayer 
tunneling resistance which otherwise should ideally be infinite at T = K [22t l2l] . 

A normal-metal electrode on top of a mesa can sometimes create a delusive impression 
of being an effective heat-leading or heat-spreading media. Not only the overall thermal 
conductance of such a thin film is low because of its small thickness, the thermal conductivity 
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FIG. 2: (Color online) The ab- and c-axis resistivities assumed in the calculations (solid lines). The 
c-axis resistivity pc marked by the thin circle has been adopted from Ref. [9] for T > 50 K. For 
the present calculations, pc has arbitrarily been extrapolated to low temperatures while imitating 
a tendency to saturation observed in experiments elsewhere [22H24| . The data represented by the 
dashed- and dash-dotted lines have been adopted from Refs. [23] and [23], respectively. 



of material in its bulk- and thin-film forms can be very different even at room temperature, 
e.g. gold [26] (see Fig. [3]). The mean-free path of heat carriers is limited by small grain 
size making thermal conductivity of polycrystalline thin films smaller than in bulk. In the 
present calculations, I use the thermal conductivity of gold calculated from an experimentally 
measured resistivity of a 100-nm thick gold thin film by using the Wiedemann-Franz law, 
k{T) = (2.44 X 10"* V2K-2) Tp{T). 

The thermal conductivity of the PMMA layer is adopted from Ref. [2^ ; the layer is 20 /xm 
thick. The thermal conductivity of other possible epoxies and glues is not very different from 
that of PMMA [28]. 

An unknown type of "silver glue" is used in the original paper [9] to attach the single 
crystals to substrates. Note however that the thermal conductance of such a glue layer is very 
likely to be smaller than that of the PMMA layer adopted in the present calculations. Indeed, 
admixture of either metallic and non-metallic particles to a polymer results in additional 
scattering of the polymer phonons at low temperature [29] thereby making k smaller. In 
the case of a metallic-particle filler, the electronic part of k is not high either, due to a poor 
electrical connectivity between the particles. Such composite polymers have high resistivity. 
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FIG. 3: (Color online) The thermal conductivity k{T) of bulk gold |25j and gold thin films of 
different thickness. k{T) of a silver thin film is given for comparison. k{T) of the thin films 
was calculated from the measured resistivity p{T) by using the Wiedemann-Franz law. Note the 
difference between the thermal conductivities of the thin films and the bulk gold even at room 
temperature [26] . 



~ 0.5 — 1.5 voTl cm in most cases [30] accounting for barely 10% of the total k. However, 
viscosity of a glue with filler can be a factor of 10 — 100 higher than that of the unfilled 
epoxy (e.g. Stycast 1265 [5U]) which results in a much thicker layer of glue between the 
single crystal and substrate. 

IV. NEWTON'S LAW OF COOLING 



A simple self-consistent analysis of I-V's in the presence of Joule heating can be made 
using Newton's law of cooling: 



T = To + a/y(/,T) 
V{I,T) = IR{T), 



(2) 
(3) 



where a and R are the thermal- and electrical resistances, respectively. R{T) here is assumed 
to be independent of I i.e. I-V's are linear (ohmic) in the absence of self-heating while 



a =const with respect to both current and temperature. All the non-linearities of the 
overall I-V's are assumed to be due to the self-heating alone. 

However, the thermal resistance a resulting from such an analysis is merely some average 
value not reflecting a real temperature distribution in the mesa. No wonder that a is found 
to be quite small for large mesas, a ~ 3.5 K/mW [31]. This value should be compared with 
a ~ 50 — 300 K/mW directly measured for different geometries of small mesas and zigzag 
stacks of IJJ's [32H31]. 

In the case of unshunted IJJ's, the I-V's should be non-linear even without Joule heating 
because of tunneling between CuO planes which are d-wave superconductors. The den- 
sity of states of a d-wave superconductor is linear in energy at small bias. The tunneling 
conductance S(V^) of an unshunted d-wave superconductor-insulator-superconductor (S-I-S) 
tunneling junction should then be parabolic in voltage. Moreover, T,{V) should become even 
steeper on approaching the gap anomaly at high bias. Note however that for large and high 
mesas the assumption of linear I-V's can be quite reasonable (see below). 

The postulate that the self-heating is the only origin of the I-V non-linearities [5S] does 
not stand up well to the following qualitative and simple reasoning (see also other arguments 
in Ref. [36]). 

The I-V curves of IJJ's consist of multiple hysteretic quasiparticle branches, each corre- 
sponding to an increasing number of IJJ's switching to the quasiparticle state with increasing 
current. Suppose that the I-V's are hnear without heating and a in Eq. |2]is big enough to 
account for the experimentally observed non-linearity of the first branch. Then, doubling of 
the dissipation power at the second branch must lead to an even stronger deviation from the 
linear behavior. In other words, the first and the second branches would never be similar in 
shape if the heating were that severe as to be the main reason for the I-V non-linearity. 

Experiments show however that at least the first 10 — 15 I-V branches can collapse 
onto a single curve if the voltage is normalized by the count number of the corresponding 
branch [37]. This confirms that the branches have the same shape, the self-heating is mod- 
erate, and a is quite small. Of course, high current does result in a significant rise of the 
mesa temperature and in a cardinal change of I-V's, especially in mesas with large number 
of junctions. The most typical is the "S"-shape of I-V's (back-bending) that has long time 
ago been recognized as resulting from both the self- heating and non-equilibrium effects [3H] ■ 

When the temperature changes, so do the thermal conductivities of all the materials 



involved in the heat transfer from the junctions. This means that a changes, too, and the 
simphcity of Eq. [2] gets broken. One should then solve the non-linear diffusion equation 
taking into account the temperature dependencies of the thermal conductivities and the 
particular geometry of a sample. This is the purpose of the present article. 

For large- and high mesas, the self-heating becomes essential already at currents much 
smaller than the superconducting critical current Ic- To simplify calculations, the I-V's 
can quite confidently be regarded as linear at such currents and in the absence of self- 
heating. Indeed, take a mesa S" = 50 x 300 ^vc? large and h = 1 /im high typical for 
THz-emission experiments. It has about 700 IJJ's in series. The sum-gap voltage should 
reach the maximum value Vgm ~ 40 V assuming the typical gap voltage Vg = 2A/e = 60 mV 
per junction. The experimental voltage range is ~ 20 times smaller than Vgm due to self- 
heating and semiconductor-like R{T). The mesa resistance R ai V ^ can be estimated 
from the c-axis resistivity pc ~ 600 i7cm, Rsg = Pch/S ~ 400 Q. Assuming that the sub-gap 
resistance does not change with voltage (i.e. the quasiparticle branch is linear), the current 
Igm corresponding to V^ = Vgm is ~ 40 V/400 Q = 100 mA. The range of currents used in the 
typical THz-emission measurements is at least two times smaller than Igm and corresponds 
to voltages ^ Vgm/2. In other words, it is enough to assume a linear J(\^)-dependence up 
to roughly half the heating-free gap voltage to justify the use of Newton's law of cooling. 
The I-V's of the small- and shallow mesas that are least affected by the self-heating indicate 
that this simplification is acceptable jSH]- 

Summarizing the discussion above, I will keep the reasonable assumption of linear I- 
V's while let the thermal resistance a vary with temperature in modified Newton's law of 
cooling. It should be noted that Newton's law of cooling is important only for calculation 
of the self-consistent I-V's. The non-uniform temperature distribution can qualitatively be 
obtained even in the simplest case of uniform power dissipation and without any assumptions 
regarding the mesa resistance. 

V. RESULTS OF SIMULATIONS AND DISCUSSION 

The Joule heating in mesas, however trivial it may look, has not been paid due attention 
to until several spectroscopic studies of the pseudogap involving IJJ's appeared |371 liO] . 
These measurements required rather high bias current to reach the I-V regions beyond 
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the superconducting gap voltage. In the case of significant Joule heating, a normal-state 
resistance that decreases with temperature can result in the appearance of maxima in the 
S(y) = dl / dV (y)-cuTYes, resembling the notorious pseudogap humps [321111]. 

In general, the high-temperature superconductors have low thermal conductivity. In 
Bi2212, the relatively high superconducting critical current density in the c-axis direction- 
(jc ~ 200 — 10000 A cm~^) and the large superconducting energy gap (Vg ~ 50 meV) set 
the scales for current and voltage. The product of the two determines the heating power 
density P in a single IJJ: P = 10 — 500 W cm~^. 

To get the total power density, P should be multiplied by the number of junctions A^ 
in the mesa. A^ is several hundred in the case of high (/i ~ 1 fim) mesas showing THz 
radiation. An efficient removal of dissipated power becomes very important. It is worth 
noting that it is not the scary-high power densities alone that determine the temperature 
inside the mesas. The mesa geometry, the material of the substrate, the glue, the contacts 
etc. are also important parameters of the problem which have been taken into account in 
the present calculations. 

Two-dimensional temperature distributions calculated for several current densities in 
mesas w = 50- and 100 yum wide are shown in Fig. Ill a) and b), respectively. The dis- 
tributions correspond to a horizontal plane at half the mesa height. The temperature is 
spatially non-uniform and increases with current. For the high current densities, the tem- 
perature can even exceed T^ in the middle part of the mesa. The thick contour lines indicate 
the boundary between the superconducting and normal-metal regions. The normal-metal 
part grows with current and the superconducting mesa eventually breaks up into two smaller 
ones connected by the common electrode and the normal-state region. The latter can provide 
an effective resistive shunting for the rest of the mesa which will be discussed below. 

Comparison of the two panels in Fig. |4] also demonstrates that the wider the mesa the 
smaller the current density at which the normal-state region first appears. However, the 
maximum temperature attained at the maximum bias is not much different in all cases. 

Fig. |5] shows a side-by-side comparison of the temperature-distribution patterns obtained 
in the present calculations and recent laser- microscopy experiments |42|. Experimental pic- 
ture shows top-electrode voltage variations due to additional heating from a small laser spot 
which is scanned over the mesa structure. The spot is about two microns in diameter and 
is estimated to raise the local temperature in the spot by a few kelvins |12]. The voltage 
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FIG. 4: (Color online) The temperature distribution in mesas 50 ^um (a) and 100 /im (b) wide at 
different current densities. The contour lines are drawn for every 5 K. T = Tc is marked by the 
thick lines. Note that a significant part of the mesas can be overheated from the bath temperature 
To = 10 K to temperatures above Tg. 



variations are less than 0.5 mV [42j which is much smaller than the mesa voltage of 1 — 2 V. 

The laser-microscopy patterns reveal distinct rings with their diameter increasing with 
bias current. The rings resemble thick contour lines of the simulations that correspond to 
T = Tc- It is fully reasonable to expect large resistance variations at local places where 
temperature is close to Tc, giving rise to relatively larger signals than from other parts of 
the mesa in the laser-microscopy experiments. 

At some high bias current the laser microscope reveals the presence of some kind of 
standing waves giving very strong response in the experiment (see the pink spots in Fig. Isb. 
The origin of the spots is not totally clear and is now being debated. This is however beyond 
the scope of the present work and will not be discussed here 
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FIG. 5: (Color online) (a) The temperature distribution in the 50 x 300 /im^ mesa at 0.5, 0.67, 
0.8, 1.0, and 1.33 MA/cm^, from bottom to top, respectively. The color scale and contour lines are 
the same as in Fig. H (b) The experimental laser-microscopy image of the 70 x 330 /xm^ mesa at 
0.61, 0.86, 1.05, 1.15, 1.50, 1.57, and 1.96 MA/cm^ of the bias current density, from bottom to top, 
respectively. The false-color scale indicates the mesa-voltage variations due to the local additional 
heating in the laser spot which is scanned over the mesa surface. The picture has been reproduced 
from the Supplementary information of the paper by H.B. Wang, et al, Phys. Rev. Lett. 102, 
017006 (2009), Ref. 02], © American Physical Society. 

Fig. |6] shows I-V's for a mesa 50 /xm wide at different bath temperatures from 10 to 
80 K with 10-K intervals. The I-V's have been calculated self-consistently, assuming the 
three-dimensional heat- and electrical-current distributions. The dashed line indicates the 
boundary at which the maximum temperature Tmax of the temperature distribution reaches 
Tc. The curves reproduce experimental I-V's qualitatively well [9]. It is also seen that the 
back-bending decreases with Tq, also in accordance with experiments. 

Fig. [7] shows calculated I-V's for different mesas at Tq = 10 K. The triple numbers along 
each line indicate the mesa width, height, and the thickness of the underlying single crystal. 
The solid dots indicate the I-V points where Tmax = Tc. 

It is the mesa sizes that largely govern the maximum temperature whereas the thickness 
of the underlying single crystal is of minor importance. Indeed, compare curve "50-1-20" 
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FIG. 6: (Color online) The current-voltage characteristics of the 50-/ini-wide mesa at different 
bath temperatures from 10 to 80 K with 10-K intervals (solid lines). The dashed line indicates a 
boundary at which the maximum of the temperature distribution in the mesa reaches Tc. 

versus "50-1-40", and "50-1-20" versus "50-2-20" or "100-1-20". Qualitatively, the heat flow 
below the mesa into a half-space of the single crystal can be modeled as having cylindrical 
or even spherical symmetry, depending on the mesa length-to-width ratio. The thermal 
conductivity of Bi2212 is anisotropic, with kab/kc ~ 10. To roughly map the heat-transfer 
problem to an equivalent isotropic media, the effective thickness of the single crystal should 
be proportionally increased making it larger than the width of the mesa {d' ~ kab/kc d = 
200 yum > w, see Fig. IT]). In the cylindrical geometry, the largest temperature gradient is 
close to the heat source. Making the heat-transfer media around the source much thicker 
than the source itself has therefore little effect on its temperature. 

From Figs. [6] and [7] it is seen that T^ax = Tc at a bias current slightly above the current Ig 
corresponding to the maximum voltage Vg of the I-V's with negative differential conductance. 
This is in agreement with the direct measurements of the mesa temperature [32]. 

The normal-state region extends by about 5 /im deep into the pedestal of the mesa at 
high current (see Fig. Isl). This indicates that the effective cooling of the mesa can only be 
achieved if the pedestal is made thinner than that size. 
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FIG. 7: The current-voltage characteristics of the 50- ^ni- wide mesa at Tq = 10 K. The triple 
numbers indicate the mesa width, the mesa height, and the thickness of the underlying single 



crystal (in fim). The dots mark the I-V points where T„ 
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Recently, THz radiation with a relatively high output power has been observed at currents 
appreciably higher than Ig [13] as compared with the original experiments where a much 
weaker radiation was measured at I < Ig [U] . Very stable high-power radiation that extends 
to four harmonics (corresponding to 2.5 THz), all having a relatively narrow line width, 
have been observed in the experiments. 

VI. SYNCHRONIZATION OF RADIATION 



The fact that the mesa temperature can exceed Tc at / > Ig where the intense THz 
radiation has been observed makes it reasonable to suggest that the heating plays an active 
role in the synchronization of IJJ's in large mesas. 

Indeed, the normal-state part of the mesa a.t I > Ig can effectively create a shunting 
resistor for the rest of the mesa which is still in the superconducting state. Even in the 
superconducting parts, the temperature is non-uniform, varying from close to the bath 
temperature up to Tc over some 100 /im (see Fig. Isj). Both the shunting resistance and an 
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FIG. 8: (Color online) The temperature distribution across the vertical middle section of the mesa 
and the underlying single crystal at j = 1.33 MA/cm^. The dash-dotted line marks the isotherm 
T = Tc- The inset schematically shows a simplified picture of the resistively shunted section of the 
mesa. The equivalent shunting circuitry is an infinite matrix of resistors representing the in- and 
out-of-plane resistivities. 

uneven critical current density resulting from such a non-uniform temperature distribution 
can presumably help synchronization of the Josephson junctions in the mesa. 

It is known that a load impedance common to all the Josephson junctions in a series 
array can stimulate their synchronization [H]. Resistive, capacitive, and inductive loads and 
some of their combinations have been considered in Ref. [SI- In general, the overdamped 
Josephson junctions (with /3c = 27r/c-R,^C/$o ~ 1) are easier to synchronize by the common 
load. It is not straightforward to assess the equivalent Pc for IJJ's with the critical current 
density and the normal-state resistance being strong functions of coordinates. 

To estimate the effective shunting resistance resulting from the normal-state parts of 
the mesa, I divide the average voltage drop across the mesa height by the total current 
flowing through the normal-state area, /„. The voltage does not change much along the 
mesa upper surface and can be assumed to be constant. I integrate the spatially varying 
current density over the normal-state area to find J„. The effective shunting resistance 
appears to quickly decrease with bias current down to 100 — 200 Q which is equivalent to 
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0.15 — 0.3 Q per junction. This value is of the order of the normal-state resistance estimated 
for a 50 X 100 /im^-large IJJ at T < T,. 

The suggestive active role of the heating-induced resistive shunting in synchronization of 
the radiation can possibly be proven experimentally. The intrinsic resistive shunting can be 
created even for relatively small currents I < Ighy intentionally destroying superconductiv- 
ity in some parts of the mesa by, say, a local Si-ion implantation [45J. The radiation is then 
expected to occur for a wider range of bias currents and temperatures than in Refs. |9l |T3] . 

It is interesting to compare the intrinsic shunting described above with traditional tech- 
niques. A resistive shunting of IJJ's by an external resistor would imply deposition of a 
thin metal film connecting the top and bottom of the mesa. Effective resistance of the 
shunt is then decided by the relatively high contact resistance Rcont between the metal 
and superconductor. Firstly, a sufficiently low Rcont can only be obtained when the metal 
thin film is deposited on the freshly cleaved surface of Bi2212. The mesa vertical sides, as 
well as the etched down parts of the single crystal are likely to be significantly degraded 
by the lithography processing and etching thereby not allowing for low contact resistances. 
Secondly, the butt-end contact area of individual IJJ electrodes, away from the mesa top 
or bottom, is very small, of the order of the electrode thickness times the mesa width, 
~ 3A X 50 ;um = 1.5 X 10"^'* m^. Even if taking the most optimistic contact resistivity that 
has been observed in experiments, 8 x 10"^ Vt cm^ [46j, the Rcont can hardly be smaller than 
about 50 kfi. This possibly explains an earlier unsuccessful attempt to synchronize IJJ's in 
Bi2212 mesas by the external shunting resistance [17] . 

Even well below J^, the temperature is distributed non-uniformly (see Fig. 111). The spatial 
variation of jc can be a reason for synchronization of radiation, as has been suggested in 
Ref. [10]. However, the radiation should then correspond to only even frequency modes 
because of a nearly symmetric temperature- and jc variations for small bias currents. To 
account for the frequencies observed in the experiments [H], a mesa with twice the width 
should have been used [TO] . 

Nonetheless, at / > J^, when the mesa is subdivided into three parts, the symmetry along 
the mesa length gets broken. Each of the superconducting parts now has jc varying from 
zero at the side facing the overheated region and up to an almost nominal jc at the opposite 
side (see Fig. |8|. A difference with Ref. [ID] is in that jc here varies along the length of 
the mesa and it is now two effective mesas having the common load and top electrode that 
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should be considered. Although it is plausible to expect the synchronization even for this 
geometry, detailed analysis of this situation is beyond the scope of this article. 

VII. SUMMARY 

Self-heating in large-area high mesas that are used in experiments on powerful THz emis- 
sion from IJJ's is analyzed by numerically solving the diffusion equation. The current-voltage 
characteristics are calculated self-consistently showing regions of negative differential con- 
ductance. Particular mesa geometry and the temperature dependencies of all the electrical- 
and thermal materials parameters are included in the calculations. The calculations reveal 
that the mesa temperature is not only significantly higher than the bath temperature, it 
is strongly non-uniform along the mesa length and width. The middle part of the mesa 
can be overheated even above T^ at high bias that is consistent with the laser-microscopy 
experiments. The normal-state regions make up an effective resistive shunting for the rest 
of the mesa. The non-uniform temperature distribution and/or the resistive shunting can 
possibly be important for synchronization of emission from the junctions in the mesa. 
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